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We report measurements of charmed hadron production from hadronic {D — > K-k) and semilep- 
tonic (/i and e) decays in 200 GeV Au+Au collisions at RHIC. Analysis of the spectra indicates that 
charmed hadrons have a different radial flow pattern from light or multi-strange hadrons. Charm 
cross sections at mid-rapidity are extracted by combining the three independent measurements, cov- 
ering the transverse momentum range that contributes to ~90% of the integrated cross section. The 
cross sections scale with number of binary collisions of the initial nucleons, a signature of charm 
production exclusively at the initial impact of colliding heavy ions. The implications for charm 
quark interaction and thermalization in the strongly interacting matter are discussed. 
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Charm quarks are a unique tool to probe the par- 
tonic matter created in relativistic heavy-ion colhsions at 
RHIC energies. Due to their large mass (~1.3 GeV/c^), 
charm quarks are predicted to lose less energy than light 
quarks by gluon radiation in the medium In con- 
trast, recent measurements of the px distributions and 
nuclear modification factors of non-photonic electrons 
(NPE) from heavy quark decays at high show a sup- 
pression level similar to light hadrons [31 ■ This observa- 
tion renews the interest in charm production and the in- 
teractions of heavy charm quarks with the hot and dense 
matter produced in nuclear collisions at RHIC. 

Measurements of charm production at low px, in par- 
ticular radial and elliptic flow, probe the QCD medium 
at thermal scales and are thus sensitive to bulk medium 
properties like density and the drag constant or viscosity. 
Model treatments for low-px charm production, such as 
energy loss by collisional dissociation Q and in-medium 
transport using a diffusion formalism (in analog to Brow- 
nian motion) and resonance cross sections Q , can be used 
to infer transport properties such as interaction cross sec- 
tions and the medium density. Calculations of charm 
transport in strongly coupling theories using AdS/CFT 
correspondence Q may allow to further determine the 
viscosity or drag coefficient of quark-gluon matter formed 
at RHIC. Ultimately, charm quark radial flow may help 
establish whether light quarks thermalize Q . 

These considerations often assume that charm quarks 
are produced only in the early stages and their produc- 
tion rate is rehably calculable by perturbative QCD 0,[1|. 
Studies of the binary collision ( A^bin from Glauber model) 
scaling of the total charm cross section from rf-|-Au to 
Au+Au collisions can be used to test these assumptions 
and determine if charm is indeed a good probe with well- 
defined initial states. The total charm production cross 
section is also an important input in models of J/ip pro- 
duction via charm quark coalescence in a Quark Gluon 
Plasma Q. 

In this paper, we present the study of charm via mea- 
surements of very low pr muons [lO|(0.17 < pr < 0.25 
GeV/c), D° Kit at low pr {pt < 2 GeV/c), and 
NPE (0.9 < _PT < 5 GcV/c). The combination of these 
three techniques provides the most complete kinematic 
coverage (~90%) to date for charm production measure- 
ments at RHIC. They also allow the study of the charmed 
hadron spectral shape in order to explore the possibility 
of charm radial flow in Au-|-Au collisions. 

The data used for these analyses were taken with the 
STAR experiment [lH during the = 200 GeV 

Au-HAu run in 2004 at RHIC. The most central 0-80% 
of the total hadronic cross section as selected using the 
uncorrected charged particle multiplicity at mid-rapidity 
(|77| < 0.5), is used for the minimum bias (minbias) mea- 



surement. A total of 13.3 million minbias triggered events 
were used for the D° reconstruction. A separate sample 
of central events was taken using an online selection of the 
0—12% most central events based on the energy deposited 
in the two Zero-Degree Calorimeters [l2| • For muons and 
NPE analysis, 15 milhon central (0-12%) and 7.8 mil- 
lion minbias triggered events are used. All measurements 
are presented as an average of particle and anti-particle 
yields at mid-rapidity {\y\ < 1 for £)° and —1 < y < 
for /x(e) limited by the TOF detector acceptance.). 

mesons were reconstructed through their decay 
D°{D") K^TT^ with a branching ratio of 3.83%. The 
analysis is identical to that used for rf+Au collisions [l3|. 
An example invariant mass distribution after subtraction 
of the combinatorial background from mixed events is 
shown in Fig.[TJa) (full circles). The yield is extracted 
by fitting a Gaussian peak plus a linear or second-order 
polynomial to describe the residual background to the 
measured distribution (red curve in Fig. [T] (a)). The 
total systematic uncertainty on yield bin-by-bin is 
^40—50%, evaluated by varying the particle identifica- 
tion conditions and yield extraction procedures. 

Muons were identified by combining the Time Of 
Flight (TOF) and ionization energy loss jdE/dx) mea- 
sured in the Time Projection Chamber (TPC) [10|, [l^, 



Il3l | . Fig. [T]^b) shows the distribution from TOF for 
tracks with 0.17 < pt < 0.21 GcV/c. Muons arc se- 
lected in the range 0.008 < < 0.014 GeVVc^- The 
distribution of the distance of closest approach (DCA) of 
the tracks to the primary vertex is used to further sepa- 
rate yu from charm decay and from pion and kaon decays. 
Fig- WLc) shows the DCA distribution of /j, after a statis- 
tical subtraction of the DCA distribution from misiden- 
tificd pions [l^ . The remaining dominant background /i 
from it/ K weak decays have different DCA distributions 
from the prompt muons, shown in Fig.[T]^c). The prompt 
/i raw yield was obtained from a fit to the p DCA distri- 
butions with the background DCA distribution and the 
primary particle DCA distribution [lo| . Other sources 
of background [p — > 77 — > 7/x"'"/j,~, Trfiu, 

etc.) are found to be negligible from simulations using 
yields from measured spectra at RHIC. The systematic 
uncertainty is dominated by the uncertainty on the pion 
contamination and was estimated using different dE/dx 
cuts for the pions. 

Electrons are identified up to pt = 5 GeV/c by dE/dx 
from TPC after applying a TOF selection |l//3 - 1| < 



0.03 [12,1131. The 5—15% bin- by-bin systematic uncer- 



tainties in inclusive electron yields are dominated by the 
uncertainties in the raw yield extraction using different 
fit functions. The dominant sources of photonic electron 
background are photon conversions and 7r°(?7) Dalitz de- 
cays. e~^e~ pairs from these background sources were 
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FIG. 1: (a) Kaon-pion pair invariant mass distribution for 
minbias Au+Au collisions, after subtraction of mixed event 
background (solid circles) and additional subtraction of a lin- 
ear residual background (open circles). A 4a signal is ob- 
served, (b) Mass squared distribution — {p/P/j)^ from 
TOF. (c) fj, DCA distributions (open circles). Histograms 
indicate the background from tt/K decays (solid curve), the 
contribution from prompt muons (dotted curve) and the sum 
(dashed curve). 
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FIG. 2: pt distributions of invariant yields for charm- 
decayed prompt fi and NPE in different centralities. Solid 
curves are power-law combined fit for D° and leptons. Dashed 
curves are blast- wave fit. The gray bands are bin-to-bin sys- 
tematic uncertainties. 

subtracted using an invariant mass technique [l^, [l^, [l^ ■ 
The photonic background reconstruction efficiency was 
calculated in a Monte Carlo simulation of 7r° and direct 
photons according to the measured 7r+ and 7r° [l3| 
spectra, and varies from 25% at low pr to 60% at high- 
est Pt in the range studied. The systematic uncertain- 
ties (7—16%) for photonic electrons arc mainly from the 



combinatorial background uncertainties and the back- 
ground reconstruction efficiency. The remaining back- 
ground from photonic decays of 77, u, p, (j) and K was de- 
termined to be ^5% from simulations and is subtracted 
in the final result. The ratio of inclusive electrons to pho- 
tonic background increases from ~ 1 to 1.45 within the 
Pt range studied. 

Fig. [2] shows the pT spectra for D'^, ^ and NPE 
for Au-|-Au minbias events. Additional centralities are 
shown for the semi-leptonic decay measurements. We 
start with a pt spectrum function for (power-law 
or blast-wave functions) and obtain spectrum function 
from semileptonic decay for NPE and /x from that. Ad- 
ditional pT-dependent factor (upper and lower bounds 
from FONLL [8|) is applied to the Icpton function to 
take into account the contribution from bottom decays. 
These functions are used in a combined fit to ob- 
tain the mid-rapidity charm yield and (pt)- The 
from the fit was with systematic errors included. The 
correlations between systematic errors on different data 
points were also taken into account in the combined 
fit [1^. The D° (pt) as calculated from power law fits 
are 0.92±0M{stat.)±0.12{sys.) GeV/c for minbias and 
0.95±0.04±0.16 GeV/c for 0-12% central Au+Au colli- 
sions. 

The nuclear modification factors {RAuAu/dAu ) M for 
/i and NPE are shown in Fig. [3l The RAuAu/dAu are 
the ratios of the pt spectra in Au-|-Au and in rf-|-Au 
collisions appropriately scaled with the number of bi- 
nary collisions. No muon measurements were carried 
out for the d-l-Au system. The d-l-Au reference is the 
decay lepton curve from the power-law combined fit 
((pt) = 1.18±0.02±0.10GeV/c,n = 11.5±0.5±1.5) with 
£)" and NPE from previously published data [3, . The 
nuclear modification factor at low pt is consistent with 1 
to within the 25% relative uncertainty of the muon mea- 
surement, and then reduces at higher pT, reaching a value 
similar to light hadrons at high pt^^ GeV/c [l6[. Our 
result is consistent with other measurements of NPE at 
high Pt . Model calculations incorporating in-medium 
charm resonances/diffusion or collisional dissociation in 
a strongly interacting medium [1, 0, 01 j which can rea- 
sonably be applied to describe the NPE spectra down to 
low Pt, arc shown in Fig. [3l 

To study whether charmed hadrons have similar ra- 
dial fiow to light hadrons, we have included curves 
for the expected nuclear modification factor from a 
blast-wave mo del, using the freeze-out parameters for 
light hadrons [2l| (BW3 in Fig. [3]) and multi-strange 
hadrons (BW2). The data and best blast- wave fit 
(BWl) show large deviations from both these curves for 
Pt > 1 GeV/c, which suggests that the charmed hadron 
freeze-out and flow are different from light hadrons. We 
scanned the parameters to a 2-dimensional Tfo, (Pt) 
space, the results show little sensitivity to freeze-out 
temperature, but disfavor large radial flow. These flnd- 
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FIG. 3: Nuclear modification factor {RauAu/cIAu) for 0—12% 
Au+Au collisions. Bin-to-bin systematic uncertainties are 
represented by the gray bands for muons and the open boxes 
for electrons. The shaded band at unity shows the com- 
mon Nbin normalization uncertainties. Model calculations 
are presented: coalescence and fragmentation [3| (double- 
dotted curves) , and collisional dissociation of heavy meson 0] 
(double-dashed curves). The curves at pr < 2 GeV/c indicate 
blast-wave model with different freeze-out parameters (T/o in 
MeV, (A)) ^ indicated in the legend. 



ings, tog ether with the observation of large charm elliptic 
flow [23| . are consistent with the recent prediction from 
hydrodynamics 2^: elliptic flow is built up at partonic 
stage, and radial flow dominantly comes from hadronic 
scattering at later stage where charm may have already 
decoupled from the system. 

A combined fit to the NPE and ^ spectra was used to 
obtain the mid-rapidity charm cross section per nucleon- 
nucleon collision {da^J^ /dy)^ shown in Fig. 21 Yields 
are an average from the blast-wave and the power-law 
functions, which are consistent within ±10%. For min- 
bias Au-I-Au collisions, the data were included in 
the fit, assuming Njjo/Ncs = 0.54±0.05 [31 as is seen 
in p+p collisions. Fig. 2] shows da^/^ /dy as a func- 
tion of iVbin for minbias c?-|-Au, minbias and central 
Au-|-Au collisions. The charm production cross sec- 
tion at mid-rapidity scales with number of binary in- 
teractions from d-|-Au [l3| to central Au-|-Au collisions. 
The quality of this scaling can be quantified by the 

(d^/dy). 



slope of dacc/dy vs A'^bin {da^^^ /dy), which is 290 
with ±15% uncorrelated uncertainty. This indicates that 
charm quarks are produced in the early stage of rela- 
tivistic heavy-ion collisions. The total cross section cr^^ 
(extrapolated to the full rapidity using the rapidity dis- 
tribution from PYTHIA by a factor of 4.7±0.7 ^) is 
1.40±0.11±0.39 mb for central and 1.29±0.12±0.36 mb 
for minbias Au+Au collisions. The central values of the 
cross sections reported by PHENIX [2^ are a factor of 
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FIG. 4: Mid-rapidity charm cross section per nucleon-nucleon 
collision as a function of Nun in ci+ Au, minbias and 0—12% 
central Au+Au collisions. The solid line indicates the average. 
FONLL prediction is shown as a band around the central 
value (thick line) [25| . 

about two smaller than ours at all measured pr- The 
difference is approximately 1.5 times the combined un- 
certainties. The FONLL calculation [2^ is the band in 
Fig. m The ratio of the minbias data over theory calcula- 
tion is A.?,±Q.h{stat.)±l.2{syst.)t^^\{theory). The upper 
theory value reproduces our result. 

In summary, we report measurements of charmed 
hadron production at mid-rapidity from analysis of 13 ^ 
Kit , ^(e) from semileptonic charm decays in minbias and 
central Au+Au collisions at RHIC. The blast-wave fits 
and the direct comparisons of the spectra suggest that 
charmed hadrons interact with and decouple from the 
system differently from the light hadrons. The dcr^^/dy 
at mid-rapidity are extracted from a combination of the 
three measurements covering the transverse momentum 
range which contributes to ~90% of the integrated cross 
section. The total cross sections are found to scale with 
the number of binary collisions. This confirms the ex- 
pected scaling of hard production processes with binary 
interactions among incoming nucleons so that charm 
quarks can be used as a calibrated probe of the early- 
stage dynamics of the system. 
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